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Abstract: We critically analyse the data available on the reduced cross-section 
in deeply inelastic ep scattering from the HI collaboration at HERA. We use 
available data on the longitudinal structure function to deduce the nature of 
dF 2 /dlrLQ^ at different for hxed values of x near x ~ 10“^. We present the 
results in a manner which effectively isolates possible higher twist effects in the 
structure function F 2 . 

PACS Nos: 12.38.-t, 13.60.Hb 

Deeply inelastic scattering (DIS) at the ep collider at HERA has provided preci¬ 
sion data on the proton structure functions over the last few years. Data on the 
slope of the structure function, F 2 {x,Q^), was presented for the hrst time last 
year Bi- This showed a surprising dip in the slope, dF 2 /d In below a: ~ 10 
(see Fig. 1), although F 2 itself continued to show a rise towards smaller x down 
to the kinematical limit. Such a dip was not anticipated or predicted by available 
parametrisations at that time |^. There have since been intensive discussions 
on this effect using leading twist as well as higher twist contributions within a 
perturbative framework. In fact, there now exist new parametrisations B a. 
which attempt to incorporate this effect, albeit in a leading twist analysis. The 
case for higher twist effects is not yet overwhelming |P] although they occur nat¬ 
urally in an operator product approach to DIS 0. In this context, it is relevant 
to ask whether there exist other data or methods which, when combined with 
existing data on dF 2 /d\n.Q^^ can clearly indicate the presence of higher twist 
terms, which have long been poorly understood in DIS. This letter addresses this 
issue. In particular, we show that the slope of the reduced DIS cross-section, a,., 
dehned in eq. (1) below, is very sensitive to the slope of F 2 . Along with available 
data on the longitudinal structure function, F^^x^Q"^), we show that the slope 
of (Tr yields information on the nature of the higher twist content of dF 2 /dhiQ^ 
and hence that of F 2 . 

Recently, the HI collaboration [§ at HERA has measured the reduced cross- 
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section, ar, 


( 1 ) 


ar{x, = 


y 


F,{x,Q^)-^Fl{x,Q^) 


which is a clear indicator of the size of the longitudinal structure function, F^^ 
compared to that of F 2 , in the measured deep inelastic ep cross-section. Here 
y is determined by = sxy with 1"+ = 1 -|- (1 — yY, where s is the square 
of the total cm energy (which is constant at HERA) and the two kinematical 
variables, and x, are as usual the momentum transfer and Bjorken scaling 
variable respectively. 

The quantity ar has been independently measured as a function of both x 
and Q^. Hence, the slope of with respect to either x or Q^, with the other 
variable kept hxed, can be determined. Furthermore, since y lies between 0 and 
1, it is clear that Fi contributes signihcantly to only at fairly large values 
of y. For example, even when y = 0.5, the quantity multiplying Fl in ar is 
only /2 = y"^/Y j^ = 0.2. Finally, Fl = 0 exactly, at leading order, due to the 
Callan-Gross relation and is non-zero only at next-to-leading order. Hence, Fl 
is suppressed (by at least one power of as) compared to F 2 . As a consequence, 
at low y, ar is essentially determined by F 2 ; conversely, any study of Fl must 
therefore be made in the large y region for maximum sensitivity. 

The data on dF 2 /d\nQ‘^ is averaged over Q^, with the average (Q^) increasing 
with increasing x in such a way that all the data corresponds essentially to low y 
^ 0.3; hence, it is relevant to ask whether the behaviour shown in Fig. 1 persists 
at all Q^. In short, what is the dependence of dF 2 /d\i\Q^ at hxed x values ? 

We analyse the HI ar [|] data with a view to extracting such a dependence. 
In order to do this, we shall separately analyse the small y and large y ar data. 
To begin with, we recognise that the dependences of Fl and dF 2 /dln.Q'^ are 
related. This is because, from eq. (1), we have 


c _ dar 
^ ~ dhiy 


dF2 £ p £ dFs 


o — dar 
^ ~ dhiy 


Q2 


= _ f p \ f dFs 

dhix ^ •’“^dhix 


( 2 ) 


Here the factors /2 = y"^/Y^ and /i = 2y‘^{2 — y)/Yl are signihcant only for large 
y. In particular, /i ~ 1 when y = 0.7; furthermore, /i ~ 2/2 over the entire y 
range. It is observed that ar is a fairly linear function of In y at all y as can be 
seen from Fig. 2, where ar has been plotted as a function of y for some selected 
X values ranging from x = 10“^ to 10“^. Hence the slope can be obtained (at 
these different x values) from straight line hts to the ar data. 


The small y data : The ^/-derivative of ar (whether at constant x or at 
constant Q'^) is insensitive to Fl or its slope when y is small. In other words. 
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Small y data 

X 


y 

dF2/d\YiQ‘^ (= dar/d\\iy\x) 

1.4 10-^ 

3.85 

0.305 

0.371 ±0.017 

2.4 10-^ 

5.24 

0.242 

0.376 ±0.010 

4.0 10“^ 

8.76 

0.243 

0.345 ± 0.005 

6.2 10-^ 

10.63 

0.190 

0.338 ±0.010 



Large y data 

X 


Fl 

dF2/d\nQ'^ (= dar/d\ny\x ± /lEj) 

0 

1 

4^ 

8.84 

0.51 ±0.3 

0.921 ±0.32 

2.4 10-^ 

15.2 

0.35 ±0.3 

0.756 ±0.32 

4.0 10“^ 

25.3 

0.33 ±0.27 

0.675 ±0.29 

6.2 10“^ 

34.7 

0.39 ±0.28 

0.758 ±0.30 


Table 1: The slope 5F2/51nQ^ evaluated using eq. 2 for different (small) x values 
for the small and large y (= 0.7) data. Note that dcTr/d\n.y\x is constant for all 
y (see Fig. 2) and directly equals dF2/dhiQ'^ at small y. The data on Fl are 
taken from the HI Collab |Q. 


the behaviour of at low y values (and hence low for a given x value) 
directly constrains the slope of F 2 . The resulting slopes, dF 2 /d\nQ^, are shown 
in comparison with those extracted differently by ZEUS in Fig. 3. Note that 
in our calculation we have obtained the slope at hxed equal to the average 
of the corresponding ZEUS data for different values of x. All the points have 
y A. 0.3; in fact the large x {x > 10“^) data have y ^ 0.15. The error bars in our 
extraction of the slope are due only to the errors arising from the straight-line 
£t to the (Jr data while the ZEUS data include both statistical and systematic 
errors. We see that there is good agreement between the two data sets, leading us 
to conclude that fiFi <C 9F2/91nQ^, so that in eq. (2) in indeed saturated by 
dF 2 /dhiQ^ at small y. The specific values of the slope, dF 2 /d\n.Q^ for certain 
X and values have been shown in Table 1 labelled as “small y" data; we shall 
need this later on in our analysis. 

The large y data : In the large y region, the Fl contribution can no longer be 
neglected. We therefore use the value of Fl determined by the HI collaboration 
0 at various x values, for hxed (corresponding to large y = 0.7) to analyse 
the large y data. Hence the large y analysis will be restricted to y = 0.7. 
Since Fl has contributions only at NLO, its slope (dF^/dlnQ^) is rather small; 
this term is further suppressed by /2 (~ 0.45 at y = 0.7). Hence, we neglect 
the contribution of the slope of Fl in what follows. We shall comment on the 
validity of this approximation later on. Then Sx in eq. (2), along with the data 
on Fl, yields a value for dF 2 /d\n.Q‘^\x at large ?/ = 0.7 within this approximation. 
Sparse data is available for Fl in the region 1 x 10“^ ^ x 6 x 10“"^; hence we 
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can extract dF 2 /d\nQ^ only at these x valnes. These valnes (which we refer to 
as “large y" data) at different bnt at the same x values as the “small y” data, 
are shown in Table 1. The “large ?/” sample obvionsly corresponds to a larger 
than the “small y” data at a given x; however, note that the average in 
the sample we have analysed increases with x. The large error bars (mnch larger 
than that of the small y data) are essentially dne to large errors in the Fi data. 

We have therefore extracted dF 2 /dlnQ‘^ as a fnnctior0 of x in the range 1- 
6 X 10“^ . These valnes are listed in Table 1. At each x valne, we have obtained 
dF 2 /d\nQ‘^ at two different valnes corresponding to small and large y data 
as discnssed above. Note that in all cases ~ 4 GeV^, which corresponds to a 
fairly stable pertnrbative regime. 


We nse the NLO GRV (1994) 0, GRV (1998) § and MRS (1998) § 
parametrisations as typical indicators of the theoretical expectation based on 
purely twist-two perturbative DGLAP [Q evolution equations. These predict 
a primarily logarithmic dependence of F 2 on (along with a small 1/Q^ piece 
from the heavy quark contributions). This implies that the slope, dF 2 /d\n.Q^, 
at hxed x is essentially flat, with a small (positive) slope due to the charm quark 
contribution. 

In Fig. 4 we show the extracted large and small y data samples as a function 
of X along with the NLO hts from the different parametrisation sets at the same 
(x, Q"^) values as the data. Since dF 2 /d\i\Q‘^ increases with the upper points 
correspond to the large y (and hence larger Q"^) sample (See Table 1). The bigger 
error bars on these points are due to the larger uncertainties in F^. The errors 
on the small y sample arise from the errors on our hts to the slope of and 
are much smaller. (The error bars on the data are very small except at the 
edges of the kinematically accessible regions; we have included them in the hts 
to the slopes but not in the error estimates of the slope). We have also plotted 
the values for dF 2 /(averaged over Q^) obtained by ZEUS at the same 
X values, at the same average as the small y data, for comparison. 

We see that all parametrisations are consistent (to within la) with the large 
y data. In the case of the small y data, which also corresponds to the smaller 
sample, the requirements are more stringent due to the smaller error bars. The 
1998 hts are in general a better ht than the 1994 one. However, it is clear that if 
the error bars decrease due to more and better data being made available, while 
the values remain near the current central ones, a pure change in normalisation 
will not suffice to ht both the data sets. This is because the central values of 
the two data sets diher by more than 100%, while the parametrisations diher 
by less than 50%. It therefore appears that the small y (or equivalently, small 
Q"^) data are suppressed relatively more than the large y data. In particular, the 


^The number of x values is restricted by Fl data. There exists more data on Fl from HI 
that has been extracted using data on Sq. This data is consistent with the existing ones 
in the region of overlap, but we shall not use them here since the data are at slightly different 
y values. 
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separation between the data at the smallest x value (= 1.4 x 10“'^) clearly indicates 
a substantial evolution of dF 2 /dlnQ^ from ~ 4 to 9 GeV^ (See Table 1). Such 
a behaviour can be attributed to higher twist effects. To elaborate further, F 2 
can be expressed in terms of a leading twist and higher twist part, as 

F2 = Clt lnQ^ + ^ , 


where the coefficients are in general functions of x and we have included only an 
additional twist-4 piece. This results in a dependence of 9^2/9 In of the 
form 


dF2 _ ^ Cht 

ainQ2 “ 


( 3 ) 


That 9 ^ 2/9 In evaluated from the parametrisations at a given x decreases 
with Q'^ is due to threshold effects from terms like occurring in the charm 

contribution. Any further suppression of c?T2/91nQ^ with decreasing then 
arises due to higher twist effects coming from the light quark sector and can 
explain the trend of the data shown in Fig. 4. Such an effect will be particularly 
visible for smaller data, such as that corresponding to the smallest x value 
in Fig. 4: this results in the small y data being suppressed more than the large 
y data. This explains the observed agreement of the large y data with existing 
twist-2 parametrisations while the small y data at the same x value disagree by 
as much as 2 ct. This disagreement disappears with increasing x for both large 
y and small y data since also increases with x. It is rather surprising that 
higher twist effects are visible at such seemingly large Q^. However, it must be 
remembered that this is the hrst time that the slope of F 2 has been so precisely 
measured. It thus appears that the Q‘^ dependence of the slope, dF 2 /d\nQ‘^, can 
provide a more sensitive test of the dependence of F 2 and hence of the elusive 
higher twist effects in deep inelastic scattering. 

Finally, we address the issue of the slope of the longitudinal structure function, 
dFL/dlnQ"^, in eq. (2), which has been neglected in this analysis. We estimate 
the size of this contribution using the GRV (1994) parametrisation. We hnd that 
this quantity does not exceed 0.1 for any of the (x,Q^) values of interest here. 
Furthermore, its contribution at small y is suppressed because of / 2 ; /2 < 0.1 for 
the small y sample, so that this contribution never exceeds a percent. For the 
large y sample, /2 is larger, /2 ~ 0.45; however, the slope of Fi at these values 
is small enough so that the term contributes less than about 5%, which is small 
compared to the size of the error bars of this data set. Hence it is reasonable and 
consistent to ignore the contribution from this term. 

In conclusion, we have analysed a limited sample of the HERA HI data [Q 
on the reduced cross-section, cr^, along with available data on the longitudinal 
structure function, |^, in order to study the dependence of dF 2 ldhiQ'^ in 
a perturbative regime. A comparison with available twist-2 NLO parametrisations 
shows indications of large dependences in the data for dF 2 /dlnQ ‘^—larger 
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than that indicated by pnrely twist-2 behavionr. This is especially trne for small 
X valnes (~ 10“^), which also correspond to a smaller (5^(~ 4 GeV^), indicating 
that snbstantial higher twist effects may be operative here. While the effect is 
clearly marked only in the hrst data point, and the data analysed is obvionsly 
limited, the trend of the data is tantalisingly similar to that expected from higher 
twist effects. We therefore nrge a detailed analysis of the cXj, data at varions (x, Q"^) 
valnes in order to shed more light on the role of higher twists in F 2 . More data 
(and improved errors on Fl) will be needed to rehne this observation. 
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Figure 1: Preliminary data on dF 2 /dlnQ^ as a function of x from ZEUS; the 
graph is taken from reference |^. The average corresponding to each x bin is 
also shown. 
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Figure 2: The data for the reduced cross-section, ar, dehned in eq. (1), and taken 
from [^, shown as a function of y along with our straight line hts to cr^ as a 
function of In y for various x values. 








Figure 3: The slope dF 2 /dlnQ"^ as obtained from cr^ data (solid circles) as de¬ 
scribed in the text (small y sample) shown as a function of x in comparison with 
ZEUS preliminary data (crosses) |]^. The of the points equals the average of 
that of the ZEUS data and is shown in Fig. 1. 
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Figure 4: The slope dF 2 /d\n (see Table 1) shown as a function of x in compari¬ 
son with standard parametrisations, as obtained from the data given in Refs, p, Q 
(see text for details). The upper (lower) circles correspond to the large (small) y 
data with the corresponding parametrisations shown as solid (dashed) lines. The 
preliminary ZEUS data (at small y) are also shown (as crosses) in the figure. 
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